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ABSTRACT 

We present results from observations of 29 BL Lacertae objects, taken with the 
Whipple Observatory 10 m Gamma- Ray Telescope between 1995 and 2000. The 
observed objects are mostly at low redshift (z < 0.2) but observations of objects 
of z up to 0.444 are also reported. Five of the objects are EGRET sources and 
two are unconfirmed TeV sources. Three of the confirmed sources of extragalactic 
TeV gamma rays were originally observed as part of this survey and have been 
reported elsewhere. No significant excesses are detected from any of the other 
objects observed, on time scales of days, months or years. We report 99.9% 
confidence level flux upper limits for the objects for each observing season. The 
flux upper limits are typically 20% of the Crab flux although, for some sources, 
limits as sensitive as 6% of the Crab flux were derived. The results are consistent 
with the synchrotron-self-Compton (SSC) model predictions considered in this 
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work. 

Subject headings: BL Lacertae objects: general — gamma rays: observations - 
galaxies: jets 



1. Introduction 



BL Lacertae objects (BL Lacs) are members of the blazar class of active galactic nuclei 
(AGN). Like all blazars, they exhibit rapid, large amplitude variability at all wavelengths, 
high optical and radio polarization, and in some cases, apparent superluminal motion and/or 
gamma-ray emission. All of these observational properties lead to t he widely held belief tha t 
blazars are AGN with jets oriented nearly along our line of sight (jUrry fc Padovanilll995l ). 
The broadband spectral energy distribution (SED) of blazars, when plotted as uF u versus 
frequency, shows a double peaked shape, with a smooth extension from radio to between 
IR and X-ray frequencies (depending on the specific blazar type), followed by a distribution 
that typically starts in the X-ray band and can peak in the gamma-ray ban d, at energies as 



high as several hundred GeV. EGRET has detecte d more than 65 blazars (IHartman et al. 
19991 ). 14 of which have been identified as BL Lacs (IDermer fc David l2000l ). 



The low energy part of the blazar SED is believed to be incoherent synchrotron radiation 
from a relativistic electron-positron plasma in the bl azar jet. The origin of the high energy 



emission is still a matter of considerable debate (e.g.. lBuckleylll998l ; lMannheimlll998l ). Lep- 
tonic models are the most popular models used to explain the observed gamma-ray emission 
from blazars. In synchrotron-self-Compton (SSC) models, the gamma rays are produced 
through inverse Compton scattering of low ener gy photons 
duce the synchrotron emission at lower energies ( Konigll 



1992 



Dermer. Sch 



ickeiser fc Mastichiadis 1992 



(Konigll 


1981 


Maraschi, Ghisellini & Celotti 


Sikora et al. 


1994; 


Bloom & Marscherl 


1996; 



Sikora fc Madejskill200ll ). In external Compton models, the dominant source of seed pho- 
tons for upscattering in the inverse Compton process are ambient photons from the cen- 
tral accretion flow, the accretion disk, the broad line region, the torus, the local infra-red 



background, or some combination of these (ISikora et a 



1994 



Dermer fc Schlickeiser 1993 



Blandford fc Levinsonl Il995l ; iGhisellini fc Madaul Il996l ; IWagner et al.l Il995l ) . Another set of 
models proposes that the gamma rays are the result of ultra-high en ergy (E > 10 19 eV) pro - 



tons producing TeV gamma rays as proton synchrotron radi ation (lAharonian et al.ll2000l ). 
in proton induced electromag netic cascades (jMannheiml Il998l ) . or by a combination of both 
processes ( jMucke et al.l 120031 ) . 



Among blazars, BL Lacs are believed to be the best candidates for TeV emission. They 
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have weak or absent optical emissio n lines, indicating that they m ay have less TeV-absorbing 
material near the emission region (IDermer fc Schlickeiserl Il994j ). Although they show rela- 
tively low luminosities compared to Flat Spectrum Radio Q uasars (FSRQs) and Optically 
Violent Variables (OVVs), their SEDs peak at higher energies ( IFossati et al.lll998l ). In blazar 
unific ation models in which electrons are assumed to be the progenitors of the gamma rays 



Ghisellini et al.lll998l ). the lower luminosity of the BL Lacs, relative to the FSRQs 



(e.g. 

implies that in BL Lacs, electrons are cooled less efficiently. This means that the electrons 
reach higher energies and that the emitted gamma rays will therefore also be of higher energy. 

Consistent with these expectations, BL Lacs are the only type of blazars detected at 
very high energies (VHE, E > 250 GeV). Curr ently, six B L Lacs are confirm ed sources of 
VHE gamma rays Mrk421 (IPunch et alJll992h. Mrk 501 (buinn et ali ll996fl. H 1426+428 
dHoran et all l2002h and 1ES 1959+6 50 iNishivama et all |200oJ); the original detection of 
the ne arby BL Lac, 1ES .2344 +514 JCatanese et all Il998h was confirmed by the HEGRA 
grou p (iTluczvkont et al.l 120031) while that of the Southern Hem isphere BL Lac, PKS 2155 - 
304 (jchadwick et al. 1999 ) was confirmed recently by H.E.S.S. ( jDjannati-Atai et al. 2003 ). 



Unco nfirmed detections have been reported for the BL Lac objects 3C66A (jNeshpor et al. 



1998h and BL Lacertae JNeshpor et allboQlh . The TeV blazars, 1ES 2344+514, H 1426+428 
and 1ES 1959+650 were originally observed at Whipple as part of the BL Lac survey de- 
scribed here. Since they were detected during these observations, the results w ere described 
in detail elsewhere ( Catanese et al. 1998 ; Horan et al. 2002 : Holder et al.ll2003l ) and are sum- 
marized here. 



Padovani fc Giommil (11995al ) introduced the terminology "Low-frequency peak BL Lacs" 
(LBLs) to describe those BL Lacs in which the lower energy SED peak occurs in the radio 
band and "High-frequency peaked BL Lacs" (HBLs) for those whos e lower energy peak 



occurs in the X-ray band. Recently, deeper BL Lacs surveys (e.g. iPerlman et al.l Il998 



Laurent-Muehleisen et al.lll999l ; ICaccianiga et al.lll999l ) have revealed evidence for the exis- 
tence of BL Lacs with properties intermediate to those in the LBL a nd HBL classes. Indeed , 
W Comae has recently been classified as an intermediate BL Lac (jTagliaferri et al.l 120001 ) . 
This suggests that, rather than being separate subclasses, LBLs and HBLs represent the 
edges of a se quence of progres sively different BL Lacs. The term "extreme blazars" was 
introduced by iGhisellinil (119991 ) to describe those BL Lacs whose first peak extends into the 
hard X-ray b a nd. S uch objects, which lie at the end of the "blazar sequence" proposed by 
Fossati et al.l (119971 ). are good candidates for TeV emission since the second peak in their 
SEDs also lies at higher energies meaning that they can be powerful at TeV energies. 

With the exception of 3C 66A and BL Lacertae, all of the claimed and confirmed TeV 
gamma-ray emitting BL Lacs are HBLs. H 1426+428, MrkSOl, 1ES 2344+514 and, to a 
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lesser extent, 1ES 1959+650, all have very hard X-ray energy spectra and fall into the class 
of extreme blazars. All of the TeV blazars but these four were listed as detections in the 



Third EGRET catalog (IHartman et al.lll999l ). The second peak in the SED of these and 
other extreme objects can only be studied in detail at TeV energies as it lies above 100 GeV, 
where EGRET and GLAST are less sensitive. Ground-based gamma-ray telescopes therefore 
offer a unique opportunity to study this class of higher-peaked blazars. Indeed, the majority 
of blazars detected by EGRET were FSRQs whose second SED peak falls in the MeV to 
GeV band. 

In order to improve our understanding of the gamma-ray emission from BL Lacs, we 
need to detect more of these objects at very high energies. Several gr oups have published 



uppe r limits on VHE emiss i on from BL Lacs in th e last several years ([Roberts et al.lll998 



999I: IChadwick et all fl99~9l : kharonian et alJ boooh . including the VERITAS collaborat 



tion 



Kerrick et aDl995h . but these efforts have been on smaller groups of objects, or not directed 



specifically at BL Lacs. In this paper we present the results of our BL Lac observing program 
from January 1995 to July 2000. The database comprises observations of 29 objects totaling 
143 hours. We present the results of searches for emission spanning time-scales of 30 minutes 
to six years. No statistically significant excess emission above the background is found and 
we discuss the implications of the non-detections. 



T hese results supersede pre liminary analyses presented in conference proceedings (ICatanese et al 



1997al ; iHoran et al.ll2000l . 120031 ). Five of the objects presented here w ere later observed more 



i ntens ively as part of a HBL observation program carried out by Ide la Calle Perez et al. 



(120031). In that survey, eight HBL s, selected from a list of TeV candidate objects derived 
by lCostamante fc Ghisellinil (120021 ). were subjected to intensive Whipple observations during 
2001 and 2002. These objects were all predicted to be VHE emitters based upon the location 
of their synchrotron peak and on their high density seed photons for the inverse Compton 
process. No evidence for TeV emission from these objects was found during this intensive 
monitoring campaign (Ide la Calle Perez et al.ll2003l ). 



2. Source List 

Table [1] lists the objects whose observations are reported here. In the table, we provide 
the object name (or names), its equatorial coordinates, its redshift and classification. The lo- 
cations of these BL Lacs, along with those of 1ES 2344+514, H 1426+428 and 1ES 1959+650, 
are plotted in Galactic co-ordinates in Figure [H 

The majority of the objects in Tabled] were selected for observations as part of three BL 
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Lac campaigns. The first was a survey of all known (circa 1995) BL Lac objects with redshift 
of < 0.1. The goal was to investigate the intrinsic characteristics of BL Lac objects that 
led to the production of TeV gamma rays. We selected low redshift objects to minimize the 
effect of the attenu ation of the VHE gamma-ray signal by pair product ion with extragalactic 
background light (jSteckerl 1 19991 ; iPrimack et al.l 1 19991 ; IVassilievI 120001 ) so that the intrinsic 
features could be compared. The second campaign was to search for TeV emission from 
HBLs in the redshift range from 0.1 to 0.2. We believed that because the HBLs have SEDs 
similar to the confirmed TeV emitters that they would be stronger TeV candidates than the 
LBLs. Because of the potential attenuation by the EBL and the larger pool of objects in the 
z=0. 1-0.2 range, w e needed to be mor e selective in our surveys. The third campaign was a 
"snapshot survey" flD'Vali et al.lll999l ) in which many BL Lacs which were considered likely 
candidates for TeV emission based on the same criteria as the second survey were observed. 
When in a flaring state, a 10 minute observation of Mrk421 or Mrk501 was enough to 
achieve a significant detection. The selected snapshot survey targets were therefore observed 
for 10 minutes each on a regular basis in the hope of catching one of them in such a flaring 
state. The objects were grouped based upon their proximity to each other on the celestial 
sphere. Objects in each group were then observed consecutively so as to minimize telescope 
slewing time. The remaining objects with known redshift were chosen because they were 
EGRET sources (W Comae, PKS 0829+046, S4 0954+65, 3C66A) or because they were a 
superluminal s ource (OQ530). RGB J172 5+118 (4U 1722+11) was observed because one 
measurement (jVeron-Cetty fc Veronl Il993l ) derived a redshift of z = 0.018, which would 
make it the closest known BL Lac object. The estimate was howev er, derived from one 
absorp tion line and many papers list its redshift as unknown (e.g., 
1995bh . 



Padovani fe Giommi 



3. Analysis Methods 
3.1. Telescope Configurations 



The VHE obse rvations reported in this paper were made with the atmospheric Cherenkov 
imaging technique ( Cawley &: Weekes 1995 ; Reynolds et al.|[l993l ) using the 10-m optical re- 
flector located at th e Whipple Observatory on Mt. Hopkins in Arizona (elevation 2.3 km) 
( jCawley et al.lll990l ). A camera, consisting of an array of photomultiplier tubes (PMTs) 
mounted in the focal plane of the reflector, records images of atmospheric Cherenkov radia- 
tion from air showers produced by gamma rays and cosmic rays. The observations reported 
here span five years and the camera of the Whipple gamma-ray telescope changed several 
times during that period. Table [2] outlines the configurations of the camera. Light concen- 
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trators are reflective cones that are mounted in front of the PMTs to improve light collection 
efficiency an d reduce albedo. Duri ng 1999, an intelligent trigger "the Pattern Selection Trig- 
ger" (PST) (IBradbury et al.lll999l ) was installed that required three adjacent tubes to record 
a signal above a certain level within a preset window. Prior to this, a trigger was declared if 
any two PMTs in the camera recorded a signal above a certain level within a preset window. 
The PST reduces the number of triggers caused by fluctuations of the night-sky background 
and thus allows the telescope to operate at lower energies. The mirror reflectivity and trigger 
settings also changed over time. Each observing season runs approximately from September 
through June. Observations are not usually carried out in July and August because the 
monsoon season, during which lightning storms strike frequently, occurs at this time. 



3.2. Gamma-ray Selection 



We characterize each Cherenkov image using a moment analysis ( iReynolds et al.lll993l ). 
The roughly elliptical shape of the image is described by the length and the width parameters, 
and its location and orientation within the telescope field of view are given by the distance 
and a parameters, respectively. We also determine the two highest signals recorded by the 
PMTs (maxl, max2) and the amount of light in the image size. In addition, we can apply 
a cut on the third moment parameter asymmetry to select gamma-ray candidates since the 
narrower tail of the image should point back toward the source location within the FOV. The 
a parameter tests whether the major axis of the image is aligned with the putative source 
location; it does not eliminate events whose major axes are parallel with the source location 
but whose image points away from it. The asymmetry parameter is not an efficient cut for 
cameras with small fields of view because the images are often truncated. 

Because of the changes in the camera discussed above in Section 13.11 the optimum cuts 
for selecting gamma rays change with time. The cuts for different camera configurations are 
listed in Table [3j They result in different sensitivities, energy ranges, and effective areas for 
each camera. This limits our ability to combine data from different observing periods into 
single upper limits. Given the variable nature of BL Lac objects however, deriving single 
upper limits for several years of observation is of dubious benefit. Instead we quote upper 
limits for each observing period. 
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3.3. Tracking Analysis 

The observations presented here were all taken in the Tracking data collection mode 
wherein only the on-source position is tracked, in runs of 28 minute duration. To estimate 
the expected background, we use those events that pass all of the gamma-ray selection 
criteria except orientation (characterized by the a parameter). We use events with values 
of a between 20° and 65° as the background region and convert the counts to an estimated 
background within the on-source region (a < 10° or 15°; see Table [3]) by multiplying the 
number of counts by a ratio determined from observations of non-source regions taken at 
other times du r ing t he observing season. This method has been described in detail by 



Catanese et al.l (119981 ). The value of the factor that converts the off-source counts to an on- 
source background estimate varies with season due to changes in the camera sensitivity and 
field of view. The estimated values for each of the observing periods are listed in Table HI 

In the case of tracking analysis, to establish the significance (S) of an excess or of a 
deficit, we use simple error propagation: 

AU - r * AW 



v/7V on + r 2 *iV- bkd + 7V2 kd *(Ar) 2 

where N on is the number of events in the on-source region (designated by the a cut in 
Table [3]), A^kd is the number of events in the background region (a = 20° —65°), and r± Ar 
is the tracking ratio and its statistical uncertainty. 



3.4. Flux Upper Limit Estimation 



After we select gamma-ray candidate events, we determine the significance of any excess 
or deficit in the observations. If the excess or the deficit is not statistically significant, as 
is the case for all observations reported here, we calc ulate a 99.9% confidence level (C.L.) 
upper limit on the count rate by using the method of iHelend (119831 ). To convert these flux 
upper limits to absolute fluxes, we first express them as a fraction of the Crab Nebula count 
rate by using observations from the same observing period. Although this method assumes 
a Crab-like spectrum for the BL Lacs, it corrects for season to season variations in factors 
like PMT gain and mirror reflectivity which affect the telescope response, and therefore its 
gamma-ray count rate. The count rates observed for the Crab Nebula for the observing 
periods reported here are given in Table HI Analysis of the Crab Nebula data shows that 
for runs taken under good we ather condit ions, the gamma-ray count rate does not change 
significantly within a season (jQuinnlll998l ). We can therefore assume that the gamma-ray 
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count rate for a source can be reliably expressed in terms of the Crab Nebula flux over the 
course of the season. 

Once we have the flux limit expressed as a fraction of the Crab Nebula count rate, we 
multiply it by the integral Crab Nebula flux (in units of photons cm~ 2 s _1 ) above the peak 
response energy of the observations (E p ). We define E p as the e nergy at which th e collection 



area folded with an E -2,5 spectrum, that of the Crab Nebula (iHillas et al.lll998l ). reaches a 
maximum. The integral fluxes from the Crab Nebula above E p for the different observations 
periods reported here are given in Table HI Upper limits are an estimate of the flux that 
could be present in the data set but not produce a significant excess. This is most accurately 
derived from the count rate because that is what determines the statistical significance of 
the excess. The Crab Nebula count rate and flux uncertainties affect only the normalization, 
so the flux upper limits quoted in terms of photons cm -2 s _1 have an uncertainty of ~25%, 
mainly from the uncertainty in the Crab Nebula photon flux. 



4. Results and Discussion 

Table [5] summarizes the results of the observations of the BL Lacs observed but not 
detected between January 1995 and July 2000 while Table [6] summarizes the results of the 
observations taken on H1426+428, 1ES1959+650 and 1ES2344+514, the three BL Lacs that 
were detected during this survey. For each target object, we list the observation exposure 
during each season, the significance of the excess or deficit of these observations, the maxi- 
mum significances for a night or month of observations and the flux upper limits expressed as 
fractions of the Crab count rate and in integral flux units (assuming a Crab-like spectrum). 
Many of the objects were observed over a number of different observing seasons resulting in 
a range of upper limits above different values of E p . No evidence for a statistically signif- 
icant excess or deficit was seen in the detected count rate from any of the objects for any 
of the time periods examined. The distribution of the significances for each object for each 
observing season is shown in Figure [H This distribution has a mean of 0.005 with a standard 
deviation of 0.976. The black curve shows the expected shape if the significances were drawn 
from a Gaussian distribution. The Kolmogorov-Smirnov test returns a 95% probability that 
the data are normally distributed. 

Table E summarizes the detections of H1426+428, 1ES1959+650 and 1ES2344+514, the 
three target BL Lacs that have been subsequently confirmed as TeV emitters. 

In order to investigate how much the flux upper limits change if the spectral index is 
not the same as that of the Crab Nebula, estimates of the flux upper limits for the BL Lacs 
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were made assuming source spectral indices of -2.2 and -2.8. The integral flux from the 
Crab Nebula above 300 GeV was used to scale the previously calculated upper limits. This 
integral flux was assumed to remain constant when the spectral index was changed. E p was 
not adjusted to account for the response of the telescope to the different input spectra but 
this effect should be very small. Table [7] lists the flux upper limits for each BL Lac for each 
observing season when these different source spectral indices were assumed. 



Costamante fc Ghisellinil (120021 ) have made predictions for the TeV flux from fourteen 



of the BL Lacs included in this paper using two different methods. In the first, an SSC 
model was used to fit the multiwavelength data gathered on each of the objects while in 
the second approach, a phenomenological description of t he average SEP of the blazars was 
derived based upon their observed bolometric luminosity (IFossati et al.lll998l ). The resulting 
flux predictions from each of the two methods were given above 300 GeV and above 1 TeV. 
In order to compare the upper limits presented here with these predictions, our upper limits 
that were not already derived above 300 GeV, were extrapolated to this energy. This was 
done by first expressing the flux upper limit as a fraction of the Crab flux at that energy, 
^BLLaci> Ep). This was then scaled to 300 GeV assuming a Crab-like spectrum. Thus, an 
upper limit on the integral flux for each BL Lac above 300 GeV for each observing season, 
FsLiac(>300GeV), was calculated: 



F B LLac{> SOOGeV) = F BLLac {> E p ) F Crab {> E } 




-1.5 



(2) 



Foab(>E p ) is the integral flux from the Crab Nebula above E p in units of photons cm" 
assuming an integral spectral index of -1.5. ^BLLac{> E p ) is the upper limit on the flux from 
the BL Lac above E„ expressed as a fraction of the Crab Nebula integral flux at this energy. 



When deriving their flux predictions, Costamante fc Ghisellini (2002J) did not take into 
account absorption of the gamma rays by the infra-red background. The predictions above 
300 GeV could therefore change by factors on the order of 5 for objects at redshifts above 
0.2. The upper limits presented here were compared with these predictions and, those 
of four BL Lacs, shown in Table [HI were found, during all seasons in which t hey were ob 



serve d, to be lower than the predicted fluxes accor ding to the Fossati approach (IFossati et al. 



19981 ) adapted in (jCostamante fc Ghisellinil 120021 ) . Those of two more BL Lacs, also listed 
in Table [91 were found to be lower during some of the observing periods. All of the up- 
per limits calculated were high er than the fluxes predicted using the one-zone, SSC model 
(jCostamante fc Ghisellinil 120021 ) . It should be noted however, that the upper limits quoted 
here pertai n only to the sp e cific p eriod during which the observations were made. As demon- 
strated by iBottcher et al.l (120021 ). spectral fitting of blazars is subject to very large uncer- 
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tainties when non-simultaneous multiwavelength data are used. Indeed, it is also shown 
that, even with the best currently available simultaneous optical - X-ray data, there is a 
very wide range in the predicted fluxes above 40 GeV. Given an observed X-ray flux, the 
predicted gamma-ray flux depends very sensitively on the model parameters and, even for 
simultaneous data, can vary by large factors due to the uncertainty in these parameters. 

In the absence of dedicated simultaneous multiwavelength data, it is difficult to use 
these data to constrain emission models. However, many of the objects surveyed here are 
monitored by the All Sky Monitor (ASM) on board the Rossi X-ray Timing Explorer. For 
many other TeV blazars, the soft X-ray flux has been seen to rise during periods when 
the gamma-ray flux was also high. For example, Mrk501 was observed to have a higher 
than average flux in the ASM during 1997, the same year during which the greatest TeV 
flaring activity detected to date was observed (jCatanese et al.l Il997bl ; iQuinn et al.lll999l ). 
In order to see if any of the BL Lacs surveyed here were particularly active in the X-ray 
band during these observations, the X-ray curves for the objects presented here that are 
monitored by ASM were analyzed. If heightened X-ray activity was detected in the absence 
of corresponding gamma-ray activity, this could have interesting consequences for emission 
models. Out of the 29 objects presented here, 25 of them are monitored on a regular basis by 
the ASM. The nightly average light curves for each of these were generated. The flux from 
each object was found to be, on average, very low (< 0.01 Crab units ) with no evidence for 
dramatic flaring or for any sustained period of X-ray activity. 

Current and future observing campaigns at the Whipple Observatory make use of BL 
Lac monitoring at X-ray wavelengths to try to predict when an object might be in a higher 
flux state and thus detectable in the VHE band. Elevated gamma-ray fluxes are often accom 



panie d by a corresponding increase in the X-ray flux (e.g. iMaraschi et al.lll999l ; iJordan et al. 



20011 ). Thus, by monitoring the X-ray activity from blazars, we can identify periods of in- 
creased activity during which the VHE flux may also be stronger. However, the relationship 
between the X-ray and gamma-ray flux has been sho wn to be complicated with gamma-ray 



flares being detected i n the absence of X- ray flares ([Holder et al 
2003h . and vice- versa jRebillot et alil2003h . 



2003; iKrawczynski et al. 



Finally, it should be noted that, although none of the objects presented here were 
found to have a statistically significant TeV flux during these observations, three of the con- 
firmed TeV blazars were detected during this survey. When 1ES2344+514, H1426+428 and 
1ES1959+650 were initially observed at Whipple, it was as part of this BL Lac campaign 
and, like the objects listed here, they were not detected. In subsequent years however, con- 
tinued monitoring with deeper exposures revealed these objects to be TeV emitters when 
in more active states, although not always detectable when in their quiescent state. Con- 



tinued VHE observ ations of the BL Lacs presen ted here, in particular those shown to have 
extreme properties (Ide la Calle Perez et al.ll2003l ). accompanied by monitoring of their X-ray 
flux level, may reveal many more of them to be TeV emitters. Indeed, since extreme BL 
Lacs, the best blazar candidates for TeV emission, have lower luminosity at all wavelengths 
than their lower energy peaked counterparts, their flux level often lies below the detection 
threshold of the current generation of imaging atmospheric Cherenkov telescopes when they 
are in quiescent state. As the next generation of ground-based gamma-ray telescopes comes 
online with their increased flux sensitivity, they will offer a unique opportunity to study this 
low luminosity class of blazar and should detect many more of these objects. 



Future, X-ray all- sky monitor experiments like LOBSTER (IBlack et al.ll2003l ) and EX- 
IST (IGrindlay et al.ll2002l ). with their improved flux sensitivity and increased bandwidth, will 
allow for more detailed monitoring of the X-ray emission from blazars thus providing valuable 
information which can be used to trigger observations at gamma-ray energies. These X-ray 
missions, coupled with the next generation of higher sensitivity VHE observatories such as 
VERITAS, HESS, MAGIC and CANGAROO, should allow both lower power extreme BL 
Lacs and lower frequency peaked blazars (LBLs and FSRQs) to be detected. 
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of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion 
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Table 1. Observed BL Lac 


Objects 












R A 




Dec. 










IN cLIIlt; 


1 J ZUUU ) 




(J2000) 




z 


Class a 


1ES 0033+595 


00 35 52 


6 


+59 50 


05 


u 


UoO 


1TRT 
1113 Li 


1ES 0145+138 


01 48 29 


7 


+14 02 


18 


n 
u 


1 0^ 
1ZO 


1TRT 


RGB J0214+517 


02 14 17 


9 


+51 44 


52 


n 
u 


U4y 


1TRT 


3C66A, lES0219+428 c,d 


02 22 39 


6 


+43 02 


08 


u 


AAA 

444 


T RT 


1ES 0229+200 


02 32 48 


4 


+20 17 


16 


u 


1 /in 

14U 


1TRT 


1H 0323+022, 1ES 0323+022 


03 26 14 





+02 25 


15 


u 


1/17 


1TRT 
1113 Li 


EXO 0706.1+5913, RGBJ0710+591 


07 10 30 





+59 08 


20 


u 


1 0^ 


11RT 
1113 L 


1ES 0806+524 


08 09 49 


1 


+52 18 


59 


u 


1 ^8 
loo 


UDT 
1113 L 


PKS 0829+046, RGB J0831+044 d 


08 31 48 


9 


+04 29 


39 


u 


i sn 

loU 


T RT 
1j131j 


1ES 0927+500 


09 30 37 


6 


+49 50 


26 


u 


1 88 
loo 


1113 L 


S4 0954+65, RGB J0958+655 d 


09 58 47 


2 


+65 33 


55 


n 
u 


OOO 


T RT 
1j131j 


1ES 1028+511 


10 31 18 


4 


+50 53 


36 


n 
u 


OOl 


TTRT 
1113 L 


1ES 1118+424 


11 20 48 





+42 12 


12 


n 
u 


1 9 /I 


TTRT 
1113 L 


Markarian 40 


11 25 36 


2 


+54 22 


57 


n 
u 


UZ 1 


TTRT 
1113 L 


Markarian 180, 1ES 1133+704 


11 36 26 


4 


+70 09 


27 


u 


D/1 K 
U4D 


TTRT 
1113 L 


1ES 1212+078 


12 15 10 


9 


+07 32 


04 


u 


loU 


TTRT 
1113 L 


ON 325, 1ES 1215+303 


12 17 52 


1 


+30 07 


01 


u 


loU 


T RT 
1j131j 


1H 1219+301, 1ES 1218+304 


12 21 21 


9 


+30 10 


37 


u 


1 89 
loZ 


TTRT 
1113 L 


W Comae, 1ES 1218+285 d 


12 21 31 


7 


+28 13 


59 


n 
u 


1 09 
1UZ 


T RT 


MS 1229.2+6430, RGB J1231+642 


12 31 31 


4 


+64 14 


18 


u 


1 vn 
1 ( U 


RRT 
11131j 


1ES 1239+069 


12 41 48 


3 


+06 36 


01 


n 
u 


1 c;n 
10U 


TTRT 
11131j 


1ES 1255+244 


12 57 31 


9 


+24 12 


40 





141 


TTrjT 

HBL 


OQ530, RGB J1419+543 


14 19 46 


6 


+54 23 


15 





151 


LBL 


4U 1722+11, RGB J1725+118 b 


17 25 04 


3 


+11 52 


15 





018 


HBL 


IZwl87, 1ES 1727+502 


17 28 18 


6 


+50 13 


10 





055 


HBL 


1ES 1741+196 


17 43 57 


8 


+ 19 35 


09 





084 


HBL 


3C371, 1ES 1807+698 


18 06 50 


6 


+69 49 


28 





051 


LBL 


BL Lacertae, 1ES 2200+420 d 


22 02 43 


3 


+42 16 


40 





069 


LBL 


1ES 2321+419 


23 23 52 


1 


+42 10 


59 





059 


HBL 



a HBL = high-frequency-peaked BL Lac object; LBL = low-frequency-peaked 
BL Lac object. 

b This object is sometimes quoted as having a redshift of 0.0 18. However, 
this is based on one absorption line (jVeron-Cetty fc Veronlll993l ) and is more 
commonly listed as having an unknown redshift. 

Unconfirmed source of TeV gamma rays. 
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Table 2. Whipple Camera Configurations 



Period 


1995/01- 
1996/12 


1997/01- 
1997/06 


1997/09- 
1998/12 


1998/12- 
1999/03 


1999/10- 
2000/07 


Number of PMTs 


109 


151 


331 


331 


379 a 


PMT spacing 


0?259 


0?259 


0?24 


0?24 


0?12 


Field of View 


3° 


3?3 


4?8 


4?8 


2?6 


Light concentrators 


yes 


yes 


no 


yes 


yes 


Pattern trigger 


no 


no 


no 


no 


yes 



a The camera consists of 379 inner tubes of FOV 0?12 diameter surrounded by 
three circular rings of PMTs (111 in all) of FOV 0?26 diameter. The outer rings of 
tubes were not used in this analysis and so, the parameters presented here pertain 
only to the inner 379 tubes. 
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Table 3. Analysis Cuts 





1995/01- 


1997/01- 


1997/09- 


1998/12- 


1999/10- 


Period 


1996/12 


1997/06 


1998/12 


1999/03 


2000/07 


maxl a 


>100 


>95 


>60 


>60 


>30 


max^ 


>80 


>45 


>40 


>40 


>30 


size* 


>400 


N.A. b 


N.A. 


N.A. 


N.A. 


length 


> 0?16 


> 0?16 


> 0?16 


> 0?16 


> 0?13 




< 0?30 


< 0?33 


< 0?50 


< 0?50 


< 0?25 


width 


> 0?073 


> 0?073 


> 0?073 


> 0?073 


> 0?05 




< 0?15 


< 0?16 


< 0?16 


< 0?16 


< 0?12 


distance 


> 0?51 


> 0?51 


> 0?51 


> 0?51 


> 0?40 




< 1?10 


< 1?17 


< 1?55 


< 1?55 


< 1?00 


alpha 


< 15° 


< 15° 


< 15° 


< 15° 


< 15° 


asymmetry 


N.A. 


> 0° 


> 0° 


> 0° 


> 0° 


length over size 


N.A. 


N.A. 


N.A. 


N.A. 


< 0.0004 



a Quantities are in units of digital counts (d.c): 1 d.c. ~ 1 photoelectron. 
b N.A. means the cut was not applied. 



-20- 



Table 4. Analysis Parameters 







Tracking 


Crab rate 


Peak Response 


Integral 


Period 


Ratio 


(7/min) 


Energy 


Crab flux a 


1995/01 - 


1995/08 


0.292 ± 0.005 


2.08 ± 0.15 


300 GeV 


1.26 


1995/10 - 


1996/07 


0.292 ± 0.004 


1.58 ± 0.05 


350 GeV 


1.05 


1996/10 - 


1996/12 


0.316 ± 0.004 


1.69 ± 0.07 


350 GeV 


1.05 


1997/01 - 


1997/06 


0.345 ± 0.005 


2.30 ± 0.10 


350 GeV 


1.05 


1998/01 - 


1998/12 


0.366 ± 0.002 


1.94 ± 0.15 


500 GeV 


0.60 


1998/12 - 


1999/03 


0.367 ± 0.004 


2.62 ± 0.26 


400 GeV 


0.84 


1999/10 - 


2000/07 


0.312 ± 0.002 


2.64 ± 0.12 


430 GeV 


0.76 



a Fluxes are quoted in units of 10 10 photons cm 2 s 1 above the corresponding 
peak response energy. 



- 21 - 



Table 5. The observation results for each object, for each observing period during which it 
was observed. For each BL Lac, the total combined significance for each period is given 
along with the maximum statistical significance seen over any one night and any month 

during that period. The flux upper limits are also presented for each observing period both 



in absolute terms (Flux units, f.u. 


, io- n 


cm 2 s" 


~ 1 ) and in 


Crab units (c.u 


•)• 






Observation 


Exp. 




Max. a 


A/To V /"V 


Flux 


Flux 


Ohiprt 


Period 


(hrs) 


a 


Month 


Nip-fit 


(c.u.) 


(f.u.) 


1ES 0033+595 




1995/12 


1.85 


-0.59 


-0.59 


0.19 


<0 


.200 


<2.10 


1ES 0145+138 


1996/10 - 


1996/11 


7 85 

1 .Ou 


-1 01 


08 

w . WO 


1.77 


<0 


.093 


^0 98 

\ W . C/O 




1998/11 - 


1998/12 


2.29 


0.22 


0.63 


0.63 


<0 


.512 


<3.50 




1998/12 - 


1999/01 


1 98 


-0 50 

W . WW 


38 

w . OO 


1.31 


<0 


.357 


<3 34 


RGB J0214+517 


1999/12 - 


2000/01 


U.U1 


99 


Sfi 


1.58 


<0 


.165 


<A 45 


3C66A 


1995/10 - 


1995/11 


8 00 

O . WW 


-2 00 

. WW 


-1 18 


0.82 


<0 


,056 


<0 59 

\ w .OU 


1ES 0229+200 


1996/11 - 


1996/12 


7 

i .00 


n 1 ^ 

U. LO 


^7 
\j.o i 


1.37 


<0 


113 


<r\ 1 Q 




1998/11 - 


1998/12 


2.30 


-1.08 


-0.78 


0.48 


<0 


,326 


<2.23 




1998/12 - 


1999/01 


1 78 

1. IO 


-0 40 


W. WW 


0.74 


<0 


,403 


\u. 1 U 


1H 0323+022 


1996/11 - 


1996/12 


10.18 


1.02 


1.02 


1.96 


<0 


,181 


<1.90 






1997/01 


91 


20 


20 

W . — w 


0.20 


<0 


.298 


<3 1 3 




1998/12 - 


1999/01 


3.18 


1.69 


1.73 


1.85 


<0 


,509 


<4.75 


EXO 0706.1+5913 




1996/12 


5 55 

o. OO 


-1 1fi 

- J. . J. W 


-1 1fi 

- J. . -LU 


0.17 


<0 


.087 


^0 91 




1997/01 - 


1997/03 


3 fi9 


7fi 


79 

W. 1 C7 


1.46 


<0 


.161 


<A fi9 






1998/11 


1.83 


-0.40 


-0.40 


1.49 


<0. 


.524 


<3.58 




1998/12 - 


1999/02 


1.90 


0.07 


1.36 


1.71 


<0. 


.459 


<4.29 


1ES 0806+524 


1996/02 - 


1996/03 


5 57 


46 


52 


1.04 


<0. 


,104 


<1 09 




2000/01 - 


2000/03 


4.16 


-0.29 


0.69 


0.69 


<1. 


293 


<11.4 


PKS 0829+046 


1995/01 - 


1995/04 


1 1 07 


1.25 


1 73 


3.14 


<0. 


,117 


<1 47 


1ES 0927+500 




1996/12 


5.08 


-1.92 


-1.92 


0.39 


<0. 


,064 


<0.67 




1997/01 - 


1997/04 


5.04 


-1.03 


0.22 


1.23 


<0. 


,076 


<0.80 


S4 0954+65 


1995/02 - 


1995/03 


3.70 


-1.09 


-0.50 


-0.11 


<0. 


096 


<1.21 


1ES 1028+511 


1998/12 - 


1999/02 


4.43 


0.57 


1.36 


2.03 


<0. 


,287 


<2.68 


1ES 1118+424 


1998/02 - 


1998/04 


7.30 


-0.25 


1.29 


1.80 


<0. 


,218 


<1.49 




1998/12 - 


1999/02 


3.60 


0.27 


1.04 


1.81 


<0. 


,310 


<2.90 




2000/01 - 


2000/05 


6.97 


-0.62 


1.29 


1.61 


<0. 


116 


<1.02 


Markarian 40 


2000/01 - 


2000/04 


10.16 


2.59 


1.66 


1.56 


<0. 


206 


<1.81 


Markarian 180 


1995/01 - 


1995/04 


5.55 


-0.10 


0.45 


1.07 


<0. 


,108 


<1.36 




1995/12 - 


1996/05 


20.46 


-0.26 


1.01 


1.70 


<0. 


,105 


<1.10 






1997/01 


0.79 


-0.17 


-0.17 


-0.17 


<0. 


303 


<3.18 
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Table 5 — Continued 





Observation 


Exp. 






Max. a 


Max. 


a 


Flux 


Flux 


Obiect 


Period 


(hrs) 


a 


Month 


Night 


(c.u.) 


(f.u.) 


1ES 1212+078 




1999/02 


1.13 


0. 


.44 


0.44 


1 


.83 


<0.778 


<7.26 




2000/01 - 


2000/05 


3.70 


1 


30 


1.52 


1 


.73 


<0.321 


<2.82 


ON 325 




1999/02 


0.97 


1 


.27 


1.27 


1 


.20 


<0.882 


<8.23 




2000/01 - 


2000/05 


5.05 


0. 


.88 


1.94 


1 


.62 


<0.215 


<1.89 


1H 1219+301 


1995/01 - 


1995/05 


2.77 


2. 


.71 


2.90 


2 


.95 


<0.226 


<2.85 




1997/02 - 


1997/06 


11.27 


0. 


.99 


1.55 


1 


.97 


<0.079 


<0.83 




1998/01 - 


1998/03 


1.38 


-1 


,96 


-1.27 


-0 


.32 


<0.356 


<2.43 




1998/12 - 


1999/02 


2.94 


-0. 


.08 


0.48 


0. 


.88 


<0.296 


<2.77 




2000/01 - 


2000/04 


3.69 


0. 


.04 


1.48 


1 


.18 


<0.191 


<1.68 


W Comae 


1995/02 - 


1995/04 


14.33 


-0 


.57 


-0.11 


1 


.14 


<0.052 


<0.66 




1996/01 - 


1996/05 


15.73 


-0 


.29 


0.24 


1 


.38 


<0.055 


<0.58 




1999/01 - 


1999/02 


4.43 


-0. 


.03 


0.58 


1 


.91 


<0.312 


<2.92 




2000/01 - 


2000/04 


4.72 


-0. 


.58 


1.01 


1 


.76 


<0.148 


<1.30 


MS 1229.2+6430 


1995/02 - 


1995/04 


1.39 


1 


.32 


1.08 


1 


.08 


<0.286 


<3.60 






1999/02 


2.04 


-0. 


.76 


-0.76 


1 


.18 


<0.446 


<4.16 




2000/01 - 


2000/05 


6.01 


0. 


.35 


1.72 


1 


.72 


<0.170 


<1.50 


1ES 1239+069 


1999/01 - 


1999/02 


1.73 





.78 


0.69 


1 


.74 


<0.616 


<6.04 




2000/01 - 


2000/05 


5.08 


0, 


.11 


1.19 


1 


.42 


<0.197 


<1.73 


1ES 1255+244 


1997/02 - 


1997/05 


5.54 


1 


.19 


1.01 


1 


.35 


<0.112 


<1.18 






1998/03 


0.46 


0. 


.13 


0.13 





.13 


<1.112 


<7.60 






1999/02 


1.73 


0. 


.15 


0.15 


1 


.10 


<0.508 


<4.75 




2000/01 - 


2000/05 


4.16 


-0 


.54 


1.30 


1 


.41 


<0.164 


<1.45 


OQ530 


1995/03 - 


1995/05 


7.39 


-0 


.73 


-0.15 





.76 


<0.058 


<0.73 


4U 1722+11 


1995/04 - 


1995/05 


2.77 


-0. 


.08 


0.28 





.70 


<0.124 


<1.56 


IZwl87 


1995/03 - 


1995/04 


2.31 


-1 


.27 


-0.07 





.63 


<0.086 


<1.08 




1996/04 - 


1996/05 


2.32 


0. 


.61 


0.85 


1 


.19 


<0.150 


<1.58 


1ES 1741+196 


1996/05 - 


1996/07 


9.23 


-1 


.02 


0.46 


2 


.07 


<0.053 


<0.56 






1998/05 


0.46 


-0 


.08 


-0.08 


-0. 


.08 


<1.168 


<7.99 


3C371 


1995/05 - 


1995/06 


13.04 


0. 


.41 


0.41 


1 


.68 


<0.190 


<1.23 


BL Lacertae 




1995/07 


4.62 


1. 


.07 


1.09 


1 


.09 


<0.109 


<1.37 
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Table 5 — Continued 





Observation 


Exp. 




Max. a 


Max. a 


Flux 


Flux 


Object 


Period 


(hrs) 


a 


Month 


Night 


(c.u.) 


(f.u.) 




1995/10 - 1995/11 


39.09 


-1.48 


-0.21 


0.85 


<0.038 


<0.40 




1998/05 - 1998/06 


0.92 


0.47 


0.71 


0.71 


<1.722 


<8.02 


1ES 2321+419 


1995/10 - 1995/11 


6.42 


-1.07 


1.50 


1.50 


<0.101 


<1.06 



Table 6. The observation results for the three detected BL Lacs that were originally 
observed as part of this survey. For each observing period during which they were 
observed, prior to detection, the total combined significance is given along with the 
maximum statistical significance seen over any one night and any month during that 
period. The flux upper limits are also presented for each observing period in absolute 

terms (Flux units, f.u., 10 -11 cm -2 s _1 ). 





Observation 


Exp. 




Max. cr 


Max. a 


Flux 


Object 


Period 


(hrs) 


a 


Month 


Night 


(f.u.) 


H 1426+428 


1995/06 - 1995/07 


3.48 


2.11 


2.11 


2.05 


<2.2 




1997/02 - 1997/06 


13.16 


1.70 


2.18 


1.62 


<0.8 




1998/04 


0.87 


1.70 


1.70 


1.99 


<6.7 


1ES 1959+650 


1995/06 


7.21 


1.02 


0.93 


1.17 


<YA 




1996/05 - 1996/07 


3.25 


0.33 


0.68 


1.10 


<1.5 




1998/07 


0.16 


-0.11 


-0.11 


-0.11 


<12.6 


1ES 2344+514 


1995/10 - 1996/01 


20.50 


5.82 


6.50 


5.80 


a 



a lES2344+514 was detected during the first observing season in which it was 
observed. See Table M for the details. 
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Table 7. The flux upper limits for each object, for each observing period during which it 
was observed when input spectra with different spectral indices (a) were assumed. The flux 
upper limits are given in absolute terms (Flux units, f.u., 1CT 11 cm~ 2 s _1 ). 









Peak 




Flux 


Flux 




Observation 


± LCD [J C11DC 


n • — 2 2 


n ■ — 2 ^ 


n ■ — 2 8 


Obiect 


Period 




ff u ) 


ff u ) 


ff u ) 


1ES 0033+595 




1995/12 


350 


<2.1476 


<2.1000 


<1.9579 


1ES 0145+138 


1996/10 - 


1996/11 


350 


<1.0022 


<0.9800 


<0.9137 




1998/11 - 


1998/12 


500 


<4.0605 


<3.5000 


<2.9886 




1998/12 - 


1999/01 


400 


<3.6323 


<3.3400 


<3.0565 


RGBJ0214+517 


1999/12 - 


2000/01 


430 


<1.6103 


<1.4500 


<1.2974 


3C66A 


1995/10 - 


1995/11 


350 


<0.6034 


<0.5900 


<0.5501 


1ES 0229+200 


1996/11 - 


1996/12 


350 


<1.2170 


<1.1900 


<1.1095 




1998/11 - 


1998/12 


500 


<2.5871 


<2.2300 


< 1.9042 




1998/12 - 


1999/01 


400 


<4.0891 


<3.7600 


<3.4408 


1H 0323+022 


1996/11 - 


1996/12 


350 


<1.9431 


<1.9000 


<1.7714 






1997/01 


350 


<3.2010 


<3.1300 


<2.9182 




1998/12 - 


1999/01 


400 


<5.1657 


<4.7500 


<4.3468 


EXO 0706.1+5913 




1996/12 


350 


<0.9306 


<0.9100 


<0.8484 




1997/01 - 


1997/03 


350 


<1.7283 


<1.6900 


<1.5756 






1998/11 


500 


<4.1533 


<3.5800 


<3.0569 




1998/12 - 


1999/02 


400 


<4.6655 


<4.2900 


<3.9258 


1ES 0806+524 


1996/02 - 


1996/03 


350 


<1.1147 


<1.0900 


<1.0162 




2000/01 - 


2000/03 


430 


<12.6599 


<1 1.4000 


<10.2005 


PKS 0829+046 


1995/01 - 


1995/04 


300 


< 1.4700 


<1.4700 


<1.4700 


1ES 0927+500 




1996/12 


350 


<0.6852 


<0.6700 


<0.6247 




1997/01 - 


1997/04 


350 


<0.8181 


<0.8000 


<0.7459 


S4 0954+65 


1995/02 - 


1995/03 


300 


<1.2100 


<1.2100 


<1.2100 


1ES 1028+511 


1998/12 - 


1999/02 


400 


<2.9146 


<2.6800 


<2.4525 


1ES 1118+424 


1998/02 - 


1998/04 


500 


<1.7286 


<1.4900 


<1.2723 




1998/12 - 


1999/02 


400 


<3.1538 


<2.9000 


<2.6538 




2000/01 - 


2000/05 


430 


<1.1327 


<1.0200 


<0.9127 


Markarian 40 


2000/01 - 


2000/04 


430 


<2.0100 


<1.8100 


<1.6196 


Markarian 180 


1995/01 - 


1995/04 


300 


< 1.3600 


< 1.3600 


<1.3600 




1995/12 - 


1996/05 


350 


<1.1249 


<1.1000 


<1.0256 
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Table 7— Continued 









Peak 


Flux 


Flux 


Flux 




Observation 


Response 


a : —2.2 


a : —2.5 


a : —2.8 


Object 


Period 


Energy 


(f.u.) 


(f.u.) 


( £ \ 

(f.u.) 






1997/01 


350 


<3 2521 


<3 1800 

W • X W W W 


<2 9648 


1ES 1212+078 




1999/02 


400 


<7.8954 


<7 2600 


<6.6437 




2000/01 - 


2000/05 


430 


<3.1317 


<2 8200 


<2 5233 


ON 325 




1999/02 


400 


<8 9503 


<8 2300 


<7.5314 




2000/01 - 


2000/05 


430 


<2 0989 


<1 8900 


<1.6911 


1H 1219+301 

_1_ J. J. _J_ — X >J \ KJ\J X 


1995/01 - 


1995/05 


300 


<2 8500 


<2 8500 


<2 8500 




1997/02 - 


1997/06 


350 

WWW 


<Q 8488 


<Q 8300 

\ w . (jo ww 


<0 7738 

\ W.I 1 u(J 




1998/01 - 


1998/03 


500 


<2.8191 


<2 4300 

Zj . T: (J W W 


<2 0750 

Zj . W 1 WW 




1998/12 - 


1999/02 


400 


<3.0124 


<2.7700 


<2 5349 




2000/01 - 


2000/04 


430 


<1.8657 


<1 6800 


<1 5032 


W Comae 


1995/02 - 


1995/04 


300 


<0 6600 


<0 6600 


<0 6600 




1996/01 - 


1996/05 


350 

WWW 


^0 5939 

\ W . WC7 WZj 


^0 5800 

\ W. WO WW 


^0 5408 

\ W . WttWO 




1999/01 - 


1999/02 


400 

TlWW 


<3 1 756 

\ W . X i WW 


<1 9200 

\Zj . JZiWU 


<1 6721 




2000/01 - 


2000/04 


430 


<-\ 44 XI 

\ J. . tit: -J ( 


<\ 3000 

\ X . OU WW 


<r1 1fi39 

\ x . lUO^ 


MS 1229 2+6430 


1995/02 - 


1995/04 


300 


<3 6000 


<3 6000 


<3 6000 






1999/02 


400 


<4.5241 


<4.1600 


<3 8069 




2000/01 - 


2000/05 


430 


<1 6658 


<1 5000 


< 1.3422 


1ES 1239+069 

XJ — < X iy | \J\J 


1999/01 - 


1999/02 


400 


<6 5686 


<6 0400 


<5 5273 




2000/01 - 


2000/05 


430 


<1.9212 


< 1.7300 


< 1.5480 


1ES 1255+244 

1 J 1 kJ XZjT_/T_/ | Z.' TI TI 


1997/02 - 


1997/05 


350 


<1 2068 


<1.1800 


<1.1001 






1998/03 


500 


<8.8171 


<7.6000 


<6.4896 






1999/02 


400 


<5.1657 


<4.7500 


<4.3468 




2000/01 - 


2000/05 


430 


<1.6103 


<1.4500 


<1.2974 


OQ530 


1995/03 - 


1995/05 


300 


<0.7300 


<0.7300 


<0.7300 


4U 1722+11 


1995/04 - 


1995/05 


300 


<1.5600 


<1.5600 


<1.5600 


IZwl87 


1995/03 - 


1995/04 


300 


<1.0800 


<1.0800 


<1.0800 




1996/04 - 


1996/05 


350 


<1.6158 


<1.5800 


<1.4731 


1ES 1741+196 


1996/05 - 


1996/07 


350 


<0.5727 


<0.5600 


<0.5221 






1998/05 


500 


<9.2695 


<7.9900 


<6.8226 
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Table 7 — Continued 







Peak 


Flux 


Flux 


Flux 




Observation 


Response 


a : -2.2 


a : -2.5 


a : -2.8 


Object 


Period 


Energy 


(f.u.) 


(f.u.) 


(f.u.) 


3C371 


1995/05 - 1995/06 


300 


<1.2300 


<1.2300 


<1.2300 


BL Lacertae 


1995/07 


300 


<1.3700 


<1.3700 


<1.3700 




1995/10 - 1995/11 


350 


<0.4091 


<0.4000 


<0.3729 




1998/05 - 1998/06 


500 


<9.3043 


<8.0200 


<6.8482 


1ES 2321+419 


1995/10 - 1995/11 


350 


< 1.0840 


< 1.0600 


<0.9883 



Table 8. Summary of the detections of the three objects, H1426+428, 1ES1959+650 and 
1ES2344+514, that were detected during this survey. The observing season during which 

they were detected, the peak response energy for that season and, where available, the 
detection integral flux in units of 10 -11 photons cm -2 s -1 ) above the corresponding peak 
response energy, along with a reference to the detection paper, are given for each object. 





Observing Season 


Peak Response 


Integral 


Detection 


Object 


of Detection 


Energy 


Flux 


Paper 


H 1426+428 


2000/10 - 2001/07 


280 GeV 


2.0 ± 0.3 


Horan et al. (2002) 


1ES 1959+650 


2001/10 - 2002/07 


600 GeV 


a 


Holder et al. f2003) 


1ES 2344+514 


1995/10 - 1996/07 


350 GeV 


1.1 ± 0.4 


Catanese et al. (1998) 



a No flux was quoted in the detection paper due to difficulties in performing a spectral 
analysis because of a decrease in telescope efficiency during the course of the observations. 
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Table 9. The 30 Upper Limits scaled to 300 GeV compared with the flux estimates of 



Costamante fc Ghisellinil (120021 ) (where available). Both the flux predictions and upper 
limits are given in abs olute flux units, f.u., 10~ n cm" 2 s _1 ). As was described in the text, 



the Upper Limits from ICostamante fc Ghisellinil (120021 ) were calculated using two different 
approaches, that of Costamante (Cos) and that of Fossati (Fos). 











F f>300 GeV) 


F (> 300 GeV) 


Observation 


Exp. 


ff u ) 


(f u) 




Obiect 


Period 


fhrs) 


Fos/Cos 


Extranolated 


1ES 0033+595 




1995/12 


1.85 


2.04 / 0.25 


< 2.64 


1ES 0145+138 


1996/10 - 


1996/11 


7.85 


-/- 


< 1.23 




1998/11 - 


1998/12 


2.29 




< 6.58 




1998/12 - 


1999/01 


1.98 




< 4.60 


RGB J0214+517 


1999/12 - 


2000/01 


6.01 


5.93 / 0.07 


< 2.14 


3C66A 


1995/10 - 


1995/11 


8.00 


0.14 / - 


< 0.74 


1ES 0229+200 


1996/11 - 


1996/12 


7.85 


0.96 / 0.31 


< 1.49 




1998/11 - 


1998/12 


2.30 




< 4.19 




1998/12 - 


1999/01 


1.78 




< 5.20 


1H 0323+022 


1996/11 - 


1996/12 


10.18 


0.84 / 0.01 


< 2.39 






1997/01 


0.91 




< 3.94 




1998/12 - 


1999/01 


3.18 




< 6.56 


EXO 0706.1+5913 




1996/12 


5.55 


-/- 


< 1.15 




1997/01 - 


1997/03 


3.69 




< 2.13 






1998/11 


1.83 




< 6.73 




1998/12 - 


1999/02 


1.90 




< 5.92 


1ES 0806+524 


1996/02 - 


1996/03 


5.57 


1.36 / - 


< 1.37 




2000/01 - 


2000/03 


4.16 




< 16.80 


PKS 0829+046 


1995/01 - 


1995/04 


11.07 


-/- 


< 1.47 


1ES 0927+500 




1996/12 


5.08 


-/- 


< 0.85 




1997/01 - 


1997/04 


5.04 




< 1.00 


S4 0954+65 


1995/02 - 


1995/03 


3.70 


-/- 


< 1.21 


1ES 1028+511 


1998/12 - 


1999/02 


4.43 


0.43 / - 


< 3.70 


1ES 1118+424 


1998/02 - 


1998/04 


7.30 


-/- 


< 2.80 




1998/12 - 


1999/02 


3.60 




< 4.00 




2000/01 - 


2000/05 


6.97 




< 1.51 


Markarian 40 


2000/01 - 


2000/04 


10.16 


-/- 


< 2.68 


Markarian 180 


1995/01 - 


1995/04 


5.55 


8.50 / 0.03 


< 1.36 




1995/12 - 


1996/05 


20.46 




< 1.39 
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Table 9 — Continued 











F (>300 GeV) 


F (> 300 GeV) 


Observation 


Exp. 


(f.U.) 


(f.U.) 




Object 


Period 


(nrs) 


.bos/Cos 


Extrapolated 






1997/01 


0.79 




< 4.00 


1ES 1212+078 




1999/02 


1.13 


/ 


< 10.03 




2000/01 - 


2000/05 


3.70 




< 4.17 


ON 325 




1999/02 


0.97 


0.16 / — 


< 11.37 




2000/01 - 


2000/05 


5.05 




< 2.79 


1H 1219+301 


1995/01 - 


1995/05 


2.77 


67 / 16 


< 2.85 




1997/02 - 


1997/06 


11.27 




< 1.04 




1998/01 - 


1998/03 


1.38 




< 4.57 




1998/12 - 


1999/02 


2.94 




< 3.82 




2000/01 - 


2000/04 


3.69 




< 2.48 


W Comae 


1995/02 - 


1995/04 


14.33 


/ 

/ 


< 0.66 




1996/01 - 


1996/05 


1 5 73 

A-'-J • 1 '.J 




< 73 




1999/01 - 


1999/02 


4 43 




<r 4 (19 




2000/01 - 


2000/04 


4 79 




< 1 Q9 


MS 1229 2+6430 


1995/02 - 


1995/04 


1.39 


/ 


< 3.60 






1999/02 


2.04 




< 5.75 




2000/01 - 


2000/05 


6.01 




< 2.21 


1ES 1239+069 


1999/01 - 


1999/02 


1.73 


/ 


< 7.94 




2000/01 - 


2000/05 


5.08 




< 2.56 


1ES 1255+244 


1997/02 - 


1997/05 


5.54 


/ 


< 1.48 






1998/03 


0.46 




< 14.28 






1999/02 


1.73 




< 6.55 




2000/01 - 


2000/05 


4.16 




< 2.13 


OQ530 


1995/03 - 


1995/05 


7.39 


-/- 


< 0.73 


4U 1722+11 


1995/04 - 


1995/05 


2.77 


12.8 / 0.015 


< 1.56 


IZwl87 


1995/03 - 


1995/04 


2.31 


5.19 / 0.07 


< 1.08 




1996/04 - 


1996/05 


2.32 




< 1.98 


1ES 1741+196 


1996/05 - 


1996/07 


9.23 


3.59 / 0.29 


< 0.70 






1998/05 


0.46 




< 15.00 
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r (>o(JU GeVJ 


r [> 300 Gevj 


Observation 


Exp. 


(f.u.) 


(f.u.) 




Object 


Period 


(hrs) 


Fos/Cos 


Extrapolated 


3C371 


1995/05 - 1995/06 


13.04 


-/- 


< 1.23 


BL Lacertae 


1995/07 


4.62 


3.32 / 0.17 


< 1.37 




1995/10 - 1995/11 


39.09 




< 0.50 




1998/05 - 1998/06 


0.92 




< 22.12 


1ES 2321+419 


1995/10 - 1995/11 


6.42 


-/- 


< 1.33 



+90 




-90 



Fig. 1. — The locations of the 32 BL Lacs originally observed as part of this survey are plotted 
in Galactic co-ordinates. The three objects which were subsequently detected (1ES2344, 
H 1426+428 and 1ES 1959+650) are labeled with diamonds while the 29 objects whose upper 
limits are reported here are labeled according to their redshift. The stars mark the four 
objects lying at z < 0.05; the circles mark six objects lying between z of 0.05 -0.1; the 
crosses mark the ten objects lying between z of 0.1 and 0.15; the plus symbols mark the five 
objects lying between z of 0.15 and 0.20; the triangles mark the four objects which lie at a 
z > 0.20. 
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Significance 



Fig. 2. — The significance of the deficit or excess in the detected count rate from each of the 
29 BL Lac objects for each season during which they were observed. This distribution has a 
mean of 0.005 and standard deviation of 0.976. The black curve shows the expected shape if 
the significances were normally distributed. This curve fits the data at the 95% confidence 
level. 



